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Decaying warm dark matter and neutrino masses 
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Neutrino masses may arise from spontaneous breaking of ungauged lepton number. Because of 
quantum gravity effects tfie associated Goldstone boson - the majoron - will pick up a mass. We 
determine tfie lifetime and mass required by cosmic microwave background observations so that the 
massive majoron provides the observed dark matter of the Universe. The majoron decaying dark 
matter (DDM) scenario fits nicely in models where neutrino masses arise a la seesaw, and may lead 
to other possible cosmological implications. 
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A long-standing challenge in particle cosmology is to 
elucidate the nature of dark matter and its origin. A keV 
weakly interacting particle could provide a sizeable frac- 
tion of the critical density Pcr = 1-88 x 10~'^^h^ g/cm^ 
and possibly play an important role in structure forma- 
tion, since the associated Jeans mass lies in the relevant 
range [ij. Although we now know from neutrino oscilla- 
tion experiments that neutrinos do have mass 0], recent 
cosmological data 0] as well as searches for distortions in 
beta [3| and double beta decay spectra 0] place a strin- 
gent limit on the absolute scale of neutrino mass that 
precludes neutrinos from being viable warm dark matter 
candidates Q and from playing a direct role in structure 
formation. 

If neutrino masses arise from the spontaneous vio- 
lation of ungauged lepton number there must exist a 
pseudoscalar gauge singlet Nambu-Goldstone boson, the 
majoron 0, This may pick up a mass from non- 
perturbative gravitational effects that explicitly break 
global symmetries @. Despite the fact that the majorons 
produced at the corresponding spontaneous L-violation 
phase will decay, mainly to neutrinos, they could still 
provide a sizeable fraction of the dark matter in the Uni- 
verse since its couplings are rather tiny. 

The decaying warm dark matter particle idea is not 
new in itself. However, since early attempts [l^l, there 
have been important observational developments which 
must be taken into account in order to assess its viability. 
Specially relevant are the recent cosmological microwave 
observations from the Wilkinson Microwave Anisotropy 
Probe (WMAP) For definiteness here we adopt the 
very popular possibility that neutrino masses arise a la 
seesaw p^ . 

In the following we consider the majoron decaying 
dark matter (DDM) idea in a modified ACDM cosmo- 
logical model in which the dark matter particle is iden- 
tified with the weakly interacting majoron J with mass 
in the keV range. The majoron is not stable but de- 
cays non radiatively with a small decay rate F. In this 
DDM scenario, the anisotropics of the cosmic microwave 



background (CMB) can be used to constrain the lifetime 
T = F^^ and the present abundance flj of the majoron; 
here we show that the cosmological constraints on DDM 
majorons not only can be fulfilled but also can easily fit 
into a comprehensive global picture for neutrino mass 
generation with spontaneous violation of lepton number. 

Although majorons could result from a phase transi- 
tion, we first consider them to be produced thermally, in 
equilibrium with photons in the early Universe. In this 
case the majoron abundance nj at the present time to 
will be. owing to entropy conservation and taking into 
account their finite lifetime: 
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where i/j is the time of majoron decoupling, and Nd 
denotes the number of quantum degrees of freedom at 
that time. If Tito) > 170 GeV, then Nd = 427/4 for 
the particle content of the standard model. On the other 
hand, in the context of a supersymmetric extension of the 
SM. there would possibly be, at sufficiently early times, 
about twice that number of degrees of freedom. More- 
over, just after decoupling, the majoron to photon ratio 
r = nj{tD)/n-y{tD) is equal to 1/2. The present density 
parameter of majorons is then 
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1.25 keV 

where we have used the standard value of Nd- 

Alternatively, if majorons were produced already out 
of equilibrium there is a range of possible models, which 
we write generically as 



mj 



-*o/t 



1.25 keV^ ■ 

where the quantity f3 parametrizes our ignorance about 
both the exact production mechanism, and the exact 
value of Nd- When /3 = 1, we recover the scenario de- 
scribed above, with r = 1/2 and Nd = 427/4. 

Clearly if the majoron is to survive as a dark mat- 
ter particle it must be long-lived, tj > t^. However, a 
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more stringent bound follows by studying the effect of a 
finite majoron lifetime on the cosmological evolution and 
in particular on the CMB anisotropy spectrum. In the 
DDM scenario, due to particle decays, the dark matter 
density is decreasing faster than in the standard cosmo- 
logical picture. This changes the time t^q of radiation- 
matter equality. This means that, for a fixed f2j, there 
will be more dark matter at early times, and the equal- 
ity will take place earlier, as illustrated in Fig. [1] The 
present amount of dark matter is ^Idm = 0.25 for both 
models; = 14Gyr in the DDM model. Other rele- 
vant parameters are ujb ~ 2.23 x 10^^ and h = 0.7. The 
time at which the blue and red lines cross is the time 
of matter-radiation equality; for fixed ^dm, it shifts to 
earlier times as the majoron lifetime decreases. 

The time of matter-radiation equality has a direct ef- 
fect on the CMB power spectrum. The gravitational po- 
tentials are decaying during the radiation dominated era; 
this means that photons will receive an energy boost after 
crossing potential wells. This so-called early integrated 
Sachs- Wolfe (EISW) effect ceases when matter comes to 
dominate the Universe, since the potential are constant 
during matter domination. The overall effect is to in- 
crease the power around the first peak as the equality 
moves to later times. 
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FIG. 1: (Color online) Evolution of abundances in the stan- 
dard (thin lines) and DDM (thick lines) universe scenario: 
blue/short dashed, red/long dashed and black/solid corre- 
spond to the matter, radiation and A components, respec- 
tively. 



On the other hand, since t > to, we expect the majoron 
decays to make the gravitational potentials vary again in 
the late stage of the cosmological evolution. This will 
induce a similar effect to the one described above, only 
affecting larger scales due to the increased horizon size. 
This late integrated Sachs- Wolfe (LISW) effect results 
then in an excess of power at small multipoles. 



Both effects can be used in principle to constrain the 
majoron lifetime and cosmological abundance. In order 
to be quantitative, we have developed a modified version 



of the CAMB code [13[, which enables us to compute the 



CMB anisotropy spectrum once the majoron lifetime and 
abundance are given in addition to the standard ACDM 
model parameters. 

Even if a kcV majoron would constitute a warm dark 
matter particle, it behaves as cold dark matter insofar 
as the calculation of its effect on the CMB spectrum 
is concerned, since CMB measurements cannot discrimi- 
nate between cold and warm dark matter 14| . The latter 
behaves differently from a cold one on scales smaller than 
its free-streaming length A/^. For a particle mass in the 
kcV range, we have Xfs ^ IMpc which corresponds in 
the CMB to a multipole £ ~ few thousands. The for- 
malism needed to account for the cosmological evolution 
of an unstable relic and of its light decay products, has 
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been developed for example in Ref. |15l . Il6l | , including the 
modifications in both background quantities and pertur- 
bation evolution. 

Two distinct mechanisms effective at very different 
times characterize the effect of DDM on the CMB. It 
is therefore convenient to choose a parametrization that 
can take advantage of this fact. In particular, the "natu- 
ral" parametrization (Oj, F) has the drawback that both 
parameters affect the time of matter-radiation equality. 
It is more convenient to define the quantity 
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where is the energy density of baryons, and tcaiiy *C 
to ^ T. As long as this condition is fulfilled, the value 
of Y does not depend on the particular choice of toariy, 
since the ratio pj / Pb is asymptotically constant at small 
times. Given that teg ^ t we can use the value of Y 
to parametrize the relative abundance of majorons at 
matter-radiation equality. In order to simplify notation 
let us also define Fig = F/(10~^^sec~^); in this way, 
Fis = 1 corresponds to a lifetime r ~ 30 Gyr. 

In the parametrization (Y, F), fixing the other param- 
eters, Y determines the time of equality, while F mainly 
affects the magnitude of the LISW effect. We show in 
Fig. [5] how the two physical effects are nicely separated 
in this parametrization. We start from a fiducial model 
with Fig = and Y = 4.7; all other parameters are 
fixed to their WMAP best-fit values. The values of Fig 
and Y are chosen in such a way to give Vijh? = 0.10, so 
that this fiducial model reproduces exactly the WMAP 
best-fit. At a larger majoron decay rate Fig = 1.2, i.e., 
F~^ ~ 27 Gyr the LISW effect makes, as expected, the 
power at small multipoles increase, while the shape of 
the spectrum around the first peak does not change, 
since the abundance of matter at early times has not 
changed. Finally, increasing Y by 20% the height of the 
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first peak decreases accordingly, wliile the largest angular 
scales (small i) are nearly unaffected. A small decrease 
in power in this region is actually observed, and can be 
explained by noticing that increasing the matter content 
we delay the onset of the A dominated era, reducing the 
A contribution to the LISW effect. 

Another advantage of using the above parametrization 
is that Y is directly related to the majoron mass through: 



Y 



^ (ki) 



(5) 
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FIG. 2; (Color online) Effect of DDM parameters on the 
CMB anisotropy spectrum. The value of the parame- 
ters are as follows. Red/solid: fiducial model (Fig, Y) = 
(0, 4.7). Green/dashed: (Fig, Y) = (1.2, 4.7). Blue/dotted: 
(Fi8, Y) = (1.2, 5.6). See text. 



We are now ready to compute the constraints that CMB 
observations put on the majoron abundance and lifetime. 
As seen from Fig. [21 even a lifetime twice as long as the 
present age of the Universe, is quite at variance with re- 
spect to the WMAP data. However one must take into 
account the fact the values of the other cosmological pa- 
rameters can be arranged in such a way as to reduce or 
even cancel the conflict with observation, i. e. degen- 
eracies may be present in parameter space. In order to 
obtain reliable constraints for the majoron mass and life- 
time, we perform a statistical analysis allowing for the 
variation of all parameters. This is better accomplished 
using a Markov chain Monte Carlo approach; we used to 



this purpose the widely known COSMOMC code 

In our modified flat (f) = 1) ACDM model, all the dark 
matter is composed of majorons. This means that no sta- 
ble cold dark matter is present The 7-dimensional 
parameter space we explore therefore includes the two 
parameters (F, F) defined above, in addition to the five 
standard parameters, namely: the baryon density flbh^, 
the dimensionlcss Hubble constant /i, the rcionization op- 
tical depth r, the amplitude As and spectral index Ug 



of the primordial density fluctuations. The cosmological 
constant density J^a depends on the values of the other 
parameters due to the flatness condition. Wc compare 
our results with the CMB anisotropics observed by the 
WMAP experiment. Once the full probability distribu- 
tion function for the seven base parameters has been ob- 
tained in this way, the probability densities for derived 
parameters, such as the majoron mass mj, is obtained. 

We show our result in Fig. [3l where we give the 68% 
and 95% confidence contours in the (mj, F) plane, for 
the case (3 = 1, i.e., thermal majoron production and 
No = 427/4. We note that these parameters arc not 
degenerate one with the other, so the respective con- 
straints are independent. Similarly, we find no degen- 
eracy between Fj and the five standard parameters. The 
marginalized 1-dimensional limits for Fj and mj are: 



Fj < 1.3 X lO-^^sec"^ 
0.12 kcV < mj < 0.17 kcV 



(6) 
(7) 



Expressed in terms of the majoron lifetime our result 
implies r > 250 Gyr, nearly a factor 20 improvement with 
respect to the naive limit t > to ~ 14 Gyr, illustrating 
the power of CMB observations in constraining particle 
physics scenarios. 
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FIG. 3: (Color online) Contours of the 68% (green/dark) and 
95% (yellow/light) confidence regions in the (Fj, mj) plane. 



Let us comment on the possibility that /3 7^ 1. From eq. 
[3l it can be seen that this amounts to the transformation 
nij — > (3m J. For example, as we have already pointed 
out, Nd can be as large as 427/2, so that the above limit 
would read 0.24 keV < mj < 0.34 keV. In general, if 
we allow for the possibility of extra degrees of freedom 
in the early Universe, we always have (3 < I and then 
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mj > 0.12 keV . If instead majorons are produced non- 
thermally, one will in general have (3 > 1. 

Let us now briefly comment on the particle physics 
model. The simplest possibility is that neutrino masses 
arise a la seesaw |12| . In the basis u, v'^ (where v denote 



ordinary neutrinos, while v'^ are the SU(2) ® U(l) singlet 
"right-handed" neutrinos) the full neutrino mass matrix 
is given as 



YJv2 



(8) 



and involves, in addition to the singlet, also a Higgs 
triplet contribution [l^ whose vacuum expectation value 
obeys a "vev seesaw" relation of the type v-iV\ ~ v\. The 
Higgs potential combines spontaneous breaking of lepton 
number and of the clectrowcak symmetry. The proper- 
ties of the seesaw majoron and its couplings follow from 
the symmetry properties of the potential and were exten- 
sively discussed in [§]. Here we assume, in addition, that 
quantum gravity effects Q produce non-renormalizable 
Planck-mass suppressed terms which explicitly break the 
global lepton number symmetry and provide the majoron 
mass, which we can not reliably compute, but we assume 
that it lies in the cosmologically interesting keV range. 

In all of such models the majoron interacts mainly with 
neutrinos, proportionally to their mass leading to 
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The limits obtained above from the WMAP data can be 
used to roughly constrain the lepton number breaking 
scale as > 3 X (lO^ GeV)^ for ~ leV. 

The massive majoron has also a sub-leading radiative 
decay mode, J 77, making our DDM scenario poten- 
tially testable through studies of the diffuse photon spec- 
trum in the far ultra violet. A more extended investiga- 
tion of these schemes will be presented elsewhere includ- 
ing other cosmological data such as the large scale struc- 
ture data from the Sloan Digital Sky Survey (SDSS) 
In contrast, we do not expect the data from upcoming 
CMB experiments like Planck to substantially improve 
our bounds on the majoron decay rate, since they mainly 
affect the large angular scales where the error bars have 
already reached the limit given by cosmic variance. We 
also note that direct detection of a keV majoron is pos- 
sible in a suitable underground experiment [2ll |. 
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